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Effect of Small Irregularities on the Constitutive 
Relations for the Ionosphere 
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Irregularities in the ionosphere which are small compared with one wavelength may 
modify the constitutive relations, and hence, may affect the refractive indices for electro- 
magnetic waves. The modifications are in some ways similar to those which would be 
introduced into the Appleton-Hartree formula by a Lorentz force. The theory is given 
first for the case when the irregularities extend only in one dimension, and it is found that 
even in a loss-free medium the refractive index now has an imaginary part which might be 
associated with loss of energy from the wave by scattering. The theory for three-dimensional 
irregularities is then discussed but is more difficult, and a method of successive approxi- 
mations is used. The results indicate that small irregularities may play an important part 
in the propagation of very-low-frequency radio waves in the ionosphere. In particular, 
they may explain why " whistlers" are observed only on comparatively rare occasions. 



1. Introduction 

In otic method of deriving the formulas of the magneto-ionic theory [l], 2 the ionized 
medium is assumed to have an electric polarization P which is related to the electric field E 
through a set of equations called the "constitutive relations." By using Maxwell's equations 
together with the constitutive relations, the Appleton-Hartree formula for the refractive index 
may be deduced. The vector P is given by: 

P = Ner (1) 

where N is the number of electrons per unit volume, e is the charge on one electron, and r is 
the average displacement of an electron from its mean position. P is thus defined by taking 
an average over a volume large enough to contain many electrons. When used in Maxwell's 
equations it is assumed to represent a vector field, continuously distributed in space, and 
approximately constant over distances which are small compared to a wavelenth. There 
would be no meaning in speaking of the value of P at a specific point in the free space between 
the electrons. The use of Maxwell's equations implies a "smoothing out" process over a 
distance which must be large compared with N~ l,z . 

There is now much evidence to show that the ionosphere is an irregular medium, so that 
the electron-number density N varies from place to place. The effect of these irregularities on 
the propagation of waves through the medium may be treated by considering the energy scat- 
tered from them, or by studying their effect on the refraction and diffraction of the waves. In 
the present paper, an alternative method is suggested for studying this problem in the special 
case when the irregularities are very small compared with one wavelength. The electric field E 
and the polarization P depend on the electron-number density N and vary from place to place 
because N varies, but it is possible to find their average values, E and P , by a "smoothing out" 
process similar to that already mentioned. It is the average values which must be used in 
Maxwell's equations; and, in order to find the refractive index, it is necessary to know how P 
depends on E . A new set of constitutive relations is therefore required giving P in terms of E , 
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and this is, in general, different from the more familiar relations which give P in terms of E. 
Hence, the refractive index must also be modified when small irregularities are present. The 
purpose of the present paper is to find these modifications in some special cases. 

2. Qualitative Description of the Method 

The ionosphere is an electrically neutral medium, and it is assumed that the number of 
heavy positive ions per unit volume is the same as the number of electrons, N. Variations in N' 
from place to place, therefore, occur equally for electrons and positive ions. They are assumed 
to be statistically stationary, and the average value of N over a very large region is JV . 

Suppose that a uniform electric field, E, is applied to the ionosphere. If E varies harmoni- 
cally in time, the electron displacements also vary harmonically with the same frequency, and 
at any one instant all the electrons are displaced the same distance, r, from their mean positions 
The positive ions are assumed to be so heavy that they do not move. Let x be distance meas- 
ured in the direction of E, and suppose that N varies with x as shown in figure 1. When the 
electrons are displaced, the effect is to introduce some space charge, which is negative where N 
is decreasing in the direction of r, and positive where it is increasing. The first effect of the 
space charge is to introduce an additional electric field, AE, whose average value is zero as long 
as the electrons are all displaced by the same amount. To a first approximation, therefore, the 
average value, E , of the electric field is just the original field E, and the average value of the 
polarization is P =i\Vr. In this case, P and E are related by the same constitutive relations 
as for a homogeneous ionosphere. This results from the assumption that the electron displace- 
ments r are all the same. 

The additional electric field AE exerts a force on the electrons which changes their dis- 
placements. This in turn changes the space charge distribution, which then affects AE. This 
is a second-order effect since AE is small, but it will be shown that it causes both E and P ta 
have new average values, E and P , which are related in a different way from the relation 
between E and P. 




NO ELECTRIC FIELD 



(0) ELECTRONS UNDISPLACED 




ELECTRIC FIELD 



(b) ELECTRONS ALL DISPLACED THE SAME DISTANCE 




ELECTRIC FIELD 



'(C) ELECTRONS DISPLACEMENTS MODIFIED 



Figure 1. When electrons are displaced in an irregular 
ionized medium, space charge appears and gives rise to 
electric forces which modify the displacements. 
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3. Comparison With the Lorentz Force 

When the magneto-ionic theory was first formulated, it was believed that the electric 
force on an elect ron in i lie ionosphere was made up of two parts. One part, eE, arose from the 
electric field E impressed on the medium and the other part, 1/3 Pe (rationalized units), arose 
from the polarization of the medium itself, and is called the Lorentz force. For many crystal- 
line dielectrics, a Lorentz force must certainly be included, but there is now evidence that in 
the ionosphere it should be omitted. The strongest evidence is from the theory of "whistlers" 
which consist of very-low-frequency waves propagated in the extensive ionized region sur- 
rounding the earth, in directions close to the lines of force of the earth's magnetic field [5]. 
The Lorentz force would so modify the constitutive relations for the ionosphere, thai these 
low-frequency waves could not be propagated. 

The modification of the constitutive relations that would be introduced by small irregulari- 
ties is similar in some ways to that which would arise from the Lorentz force. Both effects 
occur because there is a force acting on the electrons which arises from the polarization P. 
For the Lorentz force this is simply proportional to P, but for small irregularities it depends 
on both P and E. For this reason, the effect of small irregularities is more complicated than 
that of a Lorentz force, and it will be shown that the additional force on the electron may 
depend upon the state of polarization of the wave. For circularly polarized waves of very low 
frequency traveling in the direction of the earth's magnetic field (whistlers), the effect of 
irregularities may have the same direction as the Lorentz force. This affects the dispersion 
law for these waves, and if the irregularities are large enough in amplitude, there is a lower 
limit to the frequencies that can be propagated. 

4. Definitions of E and D 

For a medium which is not free space, the electric intensity E is usually defined in terms 
of a long thin cavity. The component of E in a given direction is found as follows. A long 
thin cavity is imagined to be cut in the medium parallel to the given direction. Its length 
must be so small thai the electric state of the medium does not change appreciably within it. 
For waves in a homogeneous ionosphere tins means that the length is small compared to a 
wavelength. An infinitesimal test charge 8q is placed at the center of the cavity and the 
component <5F of the force acting on it in the direction of the cavity is measured. The com- 
ponent of E in this direction is 

lim 5F/5q m 

For a medium which contains discrete ions and electrons, this definition can still be used 
provided that the cross section of the cavity is very large compared with AT"*. It is the electric 
intensity so defined wmich is to be used in Maxwell's equations. The cavity definition thus 
effects the "smoothing out" process mentioned in the introduction. A similar argument applies 
to the electric displacement D which is defined in terms of the force normal to the plane of a 
flat, plate-like cavity on a test charge in the cavity. The longest dimension of the cavity 
must be small compared to a wavelength, X and the thickness must be large compared with N~K 

Now consider an ionized medium in which there are irregularities small in size compared 
to one wavelength. Suppose that the size of the irregularities is given by a length, /, which 
may be thought of as the average distance between one point of maximum electron density 
and its nearest neighbor. Then l<^<^\. The electric intensity can now be defined in two ways. 

The "local" electric intensity E is defined in terms of a long thin cavity whose length is 
very small compared with /, and whose cross section is large compared with N~\ Its value 
depends on where the cavity is placed in relation to the maximums of N, and is affected by the 
local configuration of the space charge. It is this value of E which determines the force on an 
electron. 
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The average electric intensity E is defined in terms of a long thin cavity whose length 
is very small compared with X, and whose cross section is large compared with I. In this case 
the distance of the test charge from the walls of the cavity is very large compared with I, so 
that the value of E is unaffected by local configuration of the irregularities. It is E which 
must be used in Maxwell's equations. The cavity definition here effects the "smoothing out" 
process over a distance small compared to a wavelength but large compared with I. 

In a similar way, flat plate-like cavities may be used to define the local value D and the 
average value D of the electric displacement. The local electric polarization P is given by (1). 
Its average value is denoted by P , and is not, in general, equal to N Q e r where r is the average 
value of r. The following relations may be derived by the standard methods of electromagnetic 
theory (rationalized units) : 

D = e E + P, (2) 

Do=€ E +Po. (3) 

P and E are related by the same constitutive relations as for a homogeneous medium, but 
this is not true for P and E . 

Before attempting the general theory for three-dimensional space, the theory will be 
illustrated by discussing a simpler problem in one dimension. 

5. Irregularities in One Dimension 

Suppose that the electron density N is a function only of x, and consider the effect of 
an electric field parallel to the x direction. The effect of the earth's magnetic field is neglected. 
Our problem is to calculate the effective dielectric constant of the medium. The reader may 
find it helpful to imagine that a large sample of the medium is placed between two conducting 
plates in planes x= constant, forming a parallel plate condenser whose capacity is to be found. 
D and E are both parallel to the x direction and may be written simply as D and E. The 
electron number density N is given by 

N=N +AN, (4) 

where AN is small compared with N . Similarly the displacement r of an electron in the 
x direction from its mean position is given by 

r=r +Ar. (5) 

Further 

E=E +AE. (6) 

Each of AiV, Ar, and AE is a function of x, and has an average value of zero. 
Now div D=0 so that in this case dD/dx=0 and 

D= D = constant (7) 

for all values of x. This equation automatically takes care of the effect of space charge. The 
space charge density is p=—dPjdx. It gives rise to an additional electric field which must 
be such that div D = 0. It would clearly be incorrect to write div D=p, for this would include 
the space charge twice. 

The electron displacement r is related to E thus. 

r=yE, (8) 

where 7 is a constant which is, in general, complex. If the effect of collisions is neglected, 
y=—ejmu 2 where m is the mass of an electron and co is the angular wave frequency. If the 
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average number of collisions per unit time for an electron is v, then y=el{imvu—mu 2 ). 
Now P=Ner—NeyE and hence 

D= (e +Ney)E=D = { e + (N +AN)ey} (E +AE) . (9) 

This gives 

AE=— $-. -r^j-E^ (10) 

whose average value must be zero. 

It is now necessary to make some further assumptions about the irregularities AN. Let 
AiV=£. It is assumed that £ has a Gaussian distribution about zero, which means that the 
probability, p(£) d%, that £ lies in the range £ to £+rf£ is given by 



" (i) =vi^ '"'**'' (11) 



where 



a 2 =AiV 2 . (12) 



(A quantity with a bar over it denotes the average value of that quantity.) Then the average 
value of AE is 

aE=-5= P , ?r , t e-& 2a2 d$-E, (13) 

a^hv J -co eo+eyN +eyS 

which must be zero. The integral is evaluated in appendix A. Let 

VSa^-ftS^-Q. (14) 

Then the result of equating (13) to zero is: 



Eo 

zv 

where 



-^iT^ 2 e- Q2 -2F(Q) \, (15) 

F(0)=e-° 2 f e" 2 cfa. (16) 



If electron collisions are neglected, y is real and the integrand in (13) is then real at all points 
on the real £ axis, but has a pole on this path. The principal value of the integral is then 
purely real. If y has a small imaginary part, however, the integral has an imaginary part 
which does not tend to zero as 7(y)— >0. Thus the limiting value of the integral is not the same 
as the principal value. Now it is impossible for the collision frequency to be exactly zero so 
that clearly the limiting value is the correct one to use when collisions are neglected. 
Now 

l+N ey/e =e', (17) 

where e' is the dielectric constant of a homogeneous medium with the electron number density 
N Q . The true dielectric constant is e=D o le o E . Let 



Then (14) gives 

f=(2/3)'/W-l); a=(2|8)-»V/(«'-l), (19) 

which may be substituted in (15). 
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If electron collisions are neglected: 
and (15) gives 



e'=l-X, where X=^- 72 , 



e mar 

IJ , 2 = i =^2Xl3 l ' 2 {iw 1 ' 2 e- a2 -2F(Q)}- i . 



(20) 
(21) 



This is the formula for the square of the refractive index in the one-dimensional case when 
small irregularities are present. It replaces the formula ii 2 = l— X which would be true for a 
homogeneous medium. The quantity Xis the same as Appleton's [1] x and is proportional to the 
electron density. In figure 2 curves are plotted showing how the real and imaginary parts of n 
depend upon X for various values of ft. 

Equation (21) shows that /jl 2 has an imaginary part even when electron collisions are 
neglected. Hence, some energy would be lost from an electromagnetic wave although there 
is no physical mechanism for converting electric energy into heat. It is known that irregu- 
larities cause scattering of energy from the original wave, and this must account for the imagi- 
nary part of \x\ 

When X^>1, tt is large since (3 is small. The term, e~® 2 may then be neglected, and 
F(ti) may be replaced by the first two terms of its asymptotic expansion (eq A7, app. A). 
If powers of fi higher than the first are neglected, the expression (21) becomes 



-1- 



X 

'i+fsx/(x-iy 



(22) 



which may be compared with the expression for /x 2 in a homogeneous medium when the Lorentz 
term is included, namely: 

X 



1- 



1 + (1/3)X 



(23) 
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Figure 2. The effect of small irregularities on the variation 
of R(ju) and — I(/jl) with X. 

/8 is a measure of the strength of the irregularities. For the full curves electron 
collisions are neglected. For the dashed curves, Z=v\ 'w=0.115, and j8=0. 
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The factor p/(X—l) in (22) is analogous to the Lorentz factor, 1/3, in (23). The effect of the 
irregularities has the same sign as a Lorentz force (since X^>1) } but gets progressively smaller 
as X increases. 

6. Irregularities in Three Dimensions, General Theory 

There appears to be no direct way of extending the method of the last section to deal 
with irregularities in three dimensions. In one dimension the equation, div D = 0, has the 
simple integral D= constant, but there is no correspondingly simple result for three dimensions. 
In this and the following sections, therefore, a solution by successive approximations will be 
attempted. 

Cartesian coordinates x, y, z, are used. Let i, j, k denote unit vectors parallel to the 
x, y, and z axes respectively. Then the position of a point (x, y, z) is given by the vector 
s = ix-\-}y J rkz. The electron number density N is given by (4) where AN is now an irregular 
function of x, y, z. AA^ may be expressed as the sum of many Fourier terms thus: 

AN=J1(AN) K , (24) 

K 

(AN) K =A K cos (Ks-0*), (25) 

where 

K=ui+vj+wk. (26) 

The summation is really a triple sum over all values of u, v, w. A K and <t> K are tin 4 amplitude 
and phase, respectively, of the Fourier term. 

7. First Approximation 

As a first approximation it is assumed that the electric field E in the medium is the same 
at all points and equal to E { . Then each election is displaced from its mean position by the 
same amount r, given l>\ 

ri=T-E lf (27) 

where T is a tensor which can be found from the equation of motion of an electron. Because 
of the irregularities AN, the displacement r x results in the appearance of space charge which in 
turn gives an additional electric field AEj. The contribution AE/^ to AEi from one Fourier 
component in (24) will now be found. 

A single Fourier component (AN) K is a one-dimensional variation of AiVand may be treated 
as in section 5. The additional electric field AEi (i ° must be normal to the planes of constant 
AN for this component and so is parallel to the vector K. At this stage ri is taken as a const ant, 
so that the average value of E! must remain unaffected by the irregularities. Hence, AE/^^O. 
For components parallel to K we have e E J r P=D= constant, so that e AE+AP=0. Now 
AP=e(AN) K r 1 {K) where ri (X) is the component of ri parallel to K. Hence 

AE™ = --A K cos (Ks-fe) njS-K. (28) 



--A K cos (Ks-«^§^K. 
e K z 



Here the factor {rJL)l\K\ is r x {K \ and the factor K/\K\ is a unit vector in the direction of K. 
The total additional field AE X from the first approximation is given by: 

AE^XlAEx^. (29) 

K 

8. Second Approximation 

The effect of the appearance of the additional electric field AE T is to give the electrons an 
additional displacement, 
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Ar 1 =T-AE 1 . (30) 

Both Ar x and AEi are functions of position. Ai^ gives rise to further space charge which in 
turn results in a further contribution AE 2 to the electric field. The average value of AEi is 
zero, but the average value AE 2 of AE 2 is not in general zero, and its value must now be found. 
The resultant electron displacement is i^ + Ari, which must be substituted for r! in eq (28). 
This gives: 

AE 2 <*>=~^ A K cos (K s-fe) {(r ' + y' l) ' K} -K > (31) 

and AE 2 is given by 

AE 2 =X)AE 2 ( *>. (32) 

K 

Ari, and therefore AE 2 (iC) , contain a summation over all values of K, and so (32) contains a 

double summation. Now take the average of (32). Clearly cos (KiS — <t> Kl ) cos (K 2 s— <t> K2 ) 

has the average value zero when Kit^K 2 and the average value 1/2 when K x = K 2 . The value 
of ri may be inserted from (27). Hence: 

m * ^{TK.K ? {T El .K ? . R 

4*0 K *V 

The average electric field is now given by 

E =E 1 +AE' 2 . (34) 

Next it is necessary to find the average value P of the polarization P. From eq (1): 

P=e(N +AN)(r 1 +Ar 1 ), (35) 

and hence 

P =eN r 1 +eANAr 1 . (36) 

From (28), (29), and (30): 

A2VA ri =-- S A K cos (Ks-fe) %P TKE4 cos (Ks-fe), (37) 

*0 K h. K 



and on taking the average: 
Hence from (36) and (27) 



AiVAr 1= -£ S Al ^P T K. (38) 

Ze k IV 



P„=dV„T %-f S Ah {TE,.K}TK 
Ze k IV 



(39) 



In eq (33), (34), and (39), P and E are expressed in terms of Ei. Elimination of Ei gives the 
constitutive relations between P and E . 

9. Discussion When T Is a Scalar 

If the effect of the earth's magnetic field is neglected, the tensor T reduces to the scalar y 
used in sec. 5. Then eq (33) becomes 

^-U?^{t^- k }' <40) 

The factor [(E x • K)/K 2 ] • K is simply the component of E x in the direction of K. 

It is now assumed that the irregularities are isotropic. This means that for a given value 
of \K\ the directions of the vectors,* K, are uniformly distributed in any hemisphere and the 
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probability of a given value of A K is independent of the direction of K. Select first those 
values of K for which \K\ is in a specified range, d\K\, and A K is in a specified range, dA K , and 
perform the summation for these values first. Then the last factor in (36) is simply 1/2 E t . 
Hence (40) becomes: 

— 1 yV 



AEo= T 



t el 



K 



But 



Hence (34), (41), and (42) give: 



4 K 



R 



^{^^(H)*}- 



Equation (39) becomes: 



which reduces in a similar way to 

IyeN - 



P =E, 



y 2 e 2 N$ /ANY\ 

KnJ f' 



2 Co 

The dielectric constant e is given by 

«o« E =eoE + Po- 
It is convenient to \ise the expression (17) for «', and the abbreviation (18). Then: 

e=e7{l + ^(e'-D 2 }. 
If electron collisions are neglected and (20) is used, eq (47) becomes: 

2 _ \-X 

M ~ €_ 1+^Z 2 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



which should be compared with the expressions (21) and (22) obtained for the one-dimensional 
case. The difference is not surprising since the plane sheets of space charge formed in the one- 
dimensional problem would be expected to have a larger and different effect from an irregular 
three-dimensional distribution of charge. It is probable that (48) is unreliable when X is near 
unity, for then the wave frequency is near the plasma frequency, and it is possible for the elec- 
trons to make oscillations of large amplitude. The approximations made in the derivation 
would then be invalid. This may also explain why the expression (48) is purely real. In a more 
accurate treatment it should have an imaginary part which is associated with loss of energy 
from the wave by scattering. 

10. Discussion When T is a Tensor 

The constitutive relations will now be derived for the special case of a plane wave with its 
wave normal in the direction of the earth's magnetic field, which is taken as the z axis. The 
tensor T is then given by : 

'U 2 iYU 



T= 



1 



mo?U(U 2 -Y 2 ) 



-iYU 
,0 
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U 2 





U 2 - 



(49) 



(see for example, Banerjea [2]) where U=l—iZ, Z=v/u, v is the average number of collisions 
made per unit time by an electron, Y=o) H /a), and co h I2tt is the magnetic gyrofrequency for 
electrons. This case is of interest in the study of the propagation of whistlers [5]. 

In a homogeneous medium, a wave whose normal is parallel to the magnetic field is circu- 
larly polarized, and the electric field has no component in the direction of the wave normal. 
Now Ei is the first approximation to the electric field of the wave and is, therefore, the field of 
a wave in a homogeneous medium. Hence, it has no component parallel to the z axis, and (27) 
and (49) show that r } similarly has no z component. Hence we may write: 

r 1 =ix 1 +j 2 / 1 . (50) 

There is strong evidence that irregularities in the ionosphere are elongated in the direction of 
the earth's magnetic field [4]. To allow for this, eq (26) is replaced by 

K=ui+vj+jwk, (51) 

where / is a number less than unity. It will be assumed that the irregularities of AN are isotropic 
with respect to the variables u, v, w. The effect of the factor / is as though a region of truly 
isotropic irregularities had been stretched in the z direction by a factor 1//, and / will therefore 
be called the "elongation factor". 

The expression (33) for AE 2 must now be written down for this special case. The notation 
^ij (i, j=x, y, z) will be used for elements of T. Equation (49) shows that T X2 =T yz ~T zx = 
T zy =0, T xy = — T yx . Hence the factor (T K) • K in (33) becomes 

(T K)-K=T xx u 2 +T yy v 2 +T Z2 fw 2 . (52) 

Equations (27) and (50) show that the factor { (T E x ) • K} K in (33) is given by 

(xiu+yiv) (ui+vj+fwk) . (53) 

Now when the summation in (33) is completed, terms involving the products uv, vw, wu will 
give zero, because AiVis isotropic with respect to u, v, w. Hence (53) becomes: 

{ T Ei • K } • K=*iu 2 i+^ 2 j. (54) 

Equation (49) shows further that T xx =T yv . Hence eq (33) can now be written in full using 
(52) and (54). It gives 

^ = fi ? (u*+£+f*v?)* U^*(^+^)+7VV^ (55) 

It is shown in app. B that 



V ~* (u 2 +v 2 +fw 2 ) 2 f "* (u*+v l +fw t ) a ~ 6 *g K (u 2 +v 2 +fw 2 ) 2 4 



23 A 2 K ,„, 2 . 2 , /2 „„2N2— 2] A & C.,2J_„2 I J-2.,^2— 3 23 A K /„,2 I „,2_L -f2„„2\2— 7 L (f) P N $ (56) 



U 2 W 2 



P ¥ Al (u 2 +v 2 +f 2 w 2 )~ P ? A * {u 2 +v 2 +f 2 w 2 y- M{S) ' PN *° (57) 

where L(j), M(f) are known functions of/, and /3 is given by (18). Hence (55) becomes: 



and (34) (27) give: 



AE 2 =^ 0Nl{T„Lif)+TJii(f)}r u (58) 



^=£,+^3 pNl{ T XZ L+T ZZ M}T E,. (59) 
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Next, eq (39) must be written down for this special case. By methods similar to the above 
it can be shown thai 

S A 2 K {TEl '^ } ' TK =/3i^ (60) 

where 

2 .2 

S Ak tf + tt+fW^ ^ tt 2 + 2 , 2+/ 2 w 2=Q(/) '^S (61) 

(see app B) 

Equation (60) is now substituted in (39) and (27), (50) are used. This gives 

Po=eiVoT Ex-^2 Q/3T 2 E, (62) 

which may now be combined with (59) to give the constitutive relation. 
The relations (59) and (62) are equivalent to: 

E,= ^T-'+^jNl(T xx L+T z M)\TE l , (63) 

P =eN (l-^° QP TJT E l9 (64) 

from which the vector T Ei is to be eliminated. Tins then gives 

E =RPo (65) 

where R is to be found. Now is a small quantity and its square and higher powers will be 
neglected. Then it can be shown that 

and it follows from the expression (49) for T thai 

R xz = R zx = Ryz = R zy = I 

r (67) 

The standard methods of magneto-ionic theory may now be used to find the polarization and 
refractive indices of the waves which can be propagated. The two polarizations are given 
by p=±i, so that the two waves are circularly polarized with opposite senses, and the two 
refractive indices are given by 

M 2 =l+ (R \. p , > (68) 

which leads to 

1 1 . Application to Whistlers 

Whistlers are naturally occurring radio signals propagated in the extensive ionized region 
surrounding the earth, with the wave normals in directions close to the lines of force of the earth's 
magnetic field [5]. The frequency is very low (often in the audible range) so that X and Y 
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are very large compared to unity. The waves are not heavily attenuated so that the effect 
of electron collisions must be small. Hence in (69) we take [7=1, X^>1, Y^>\. More- 
over X/Y is inversely proportional to frequency and is very large at low frequencies, so that 
X/Y^>1 [5]. One of the values of \i given t^ (69) is almost purely imaginary; the other value, 
with — Y in the denominator, refers to the wave which is believed to explain whistlers, and it 
is given by 

"' ~ —8 / X L \ 1 (70) 

Y-%X>(M-±)+±(lQX 

If the elongation factor /is not near zero, then M(f) and L(f) are both non-zero, but L/Y^^M 
since Fis large. Zmay be neglected in comparison with X 2 , so that (70) becomes: 

^ 2 ~Y-y 2 pMX 2 ' (71) 

The additional term ) f ^MX 2 in the denominator has an effect similar to that which the Lorentz 
force would have. Waves could not be propagated if 

\MX>>Y, (72) 

for then the refractive index would be purely imaginary. Hence, irregularities could modify 
the propagation of whistlers. Equation (72) shows that there is a lower limit to the frequencies 
which could be propagated, given by 

Q4PMM*y (73) 

where f N is the plasma frequency and f H is the gyrofrequency. 

If, however, the irregularities are very elongated in the direction of the earth's magnetic 
field, so that the elongation factor /is very small, then M(f) is small (see app. B), and its effect 
could be opposed by the term L/Y 2 in the denominator of (70). If / is so small that L — MY 2 , 
then the effect of irregularities on the propagation of whistlers would be small. 

It is well known that whistlers are not observed regularly, but only during comparatively 
short periods of time. They are believed to start from lightning flashes, so that sources of 
whistlers are always plentiful. The absence of whistlers for long periods must therefore be 
explained by some feature of the propagation. It is now tentatively suggested that whistlers 
cannot be observed when the path they would travel contains irregularities which are not 
greatly elongated. On the occasions when whistlers can be observed, the irregularities are 
either absent or greatly elongated in the direction of the earth's magnetic field. 

12. Conclusion 

It has been shown that irregularities in the ionosphere, which are small compared with 
one wavelength, may modify the constitutive relations and hence may affect the refractive 
indices for electromagnetic waves. The modifications are in some ways similar to those which 
would be introduced into the Appleton-Hartree formula by a Lorentz force. When the irregu- 
larities extend only in one dimension the theory can be handled by the statistical method of 
section 5, and it was found that even in a loss-free medium the refractive index now has an 
imaginary part which might be associated with loss of energy from the wave by scattering. 
For the general three-dimensional case, the theory is more difficult and has only been handled 
by a method of successive approximations. A theory similar to that of section 5 is needed for 
this case but has not so far been found. The method of section 6 onwards is less satisfactory 
because of the numerous approximations. The results, however, do indicate that small irregu- 
larities may play an important part in the propagation of very-low-frequency radio waves in 
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the ionosphere. In particular, they may explain why whistlers are observed only on compar- 
atively rare occasions. 



The author is greatly indebted to the University of Colorado, and to the Directors and 
Staffs of the High Altitude Observatory and NBS for providing him with the opportunity to 
work in Boulder. 

13. Appendix A: Evaluation of the Integral in Equation (13) 
With the notation (14), the integral in (13) becomes 

Let ^2a=w. Then (Al) gives 

Vxf J-«, w+12 



_«i2 



%L£* 



Now 



i-^r e- <r( " +0) rfT, (A3) 

i M Jo 



where the top limit is +<» when / (12) is negative, and — °° when I (12) is positive. Here / (12) 
is always positive. Hence 

£ 5Ta '-"-'J." "•""" [[_>•"-""«»] * 

-iV?r _ "e H " , '^" w ((r=2i>" (! -»'r"* e -"** 

Jo Jz8 

= -ir 1/2 [w l ' 2 e- Q2 -2F(Q)] (A4) 

where 

F(0)=e" 0, fV 2 -du. (A5) 

When (A5) is subsitituted in (A2) it leads to the result (15). 

The function F (12) has been studied by Mitchell and Zemansky [3] who give the following 
formulas : 
when 12 is small ; 

l-2^ = l_^+^-<§!+ • • • (A6) 

and when 12 is large ; 

F(Sl) ^212 L 1+ l2i2^) + l212r + (212^+ ' ' J (A7) 

14. Appendix B: Evaluation of Coefficents in Equations (55) and (60) 
According to eq (51) the vector K=iu+i*'+k/w. Its direction cosines are therefore 

7= u ™ = ____JL__ ->= f w I run 

(u 2 +v 2 +fw 2 ) 1 ^ (u 2 +v 2 +fw 2 ) 1 ' 2 ' (u 2 +v 2 +f 2 w 2 y< 2 { J 

These may be written in spherical polar coordinates: 

Z=sin 6 cos </>, ra=sin 6 sin <£, ?i=cos B. (B2) 
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In equation (55) the coefficient u A /(u 2J rv 2 +J 2 w 2 ) 2 is Z 4 =sin 4 cos 4 </>. It is assumed that 
the distribution of the directions of K is the same for all values of A K , which means that 



S -4iZ 4 =sin 4 e cos 4 S A 2 K . (B3) 



Since all values of <j> are equally probable this is equal to 



JZ ^im 4 =sin 4 6 sin 4 </> X) ^£- (B4) 

K K 



It can be shown that sin 4 </>=cos 4 <£=3/8; cos 2 c/> sin 2 <£=l/8; cos 2 </>=sin 2 </>=l/2. Hence: 



1ZAW = J^A 2 K m*=3/8 sin 4 e J2A 2 K} (B5) 



J2A 2 K l 2 rn 2 = 1/8 sin 4 6 J}A 2 K , (B6) 



S^Z W =Z)^w V=l/2 cos 2 sin 2 d J2A 2 K , (B7) 



Now let 
so that 



2>1^ 2 = ZMim 2 =l/2 sm 2 9^2 A 2 K , (B8) 

S^in 2 = cos 2 02>l| > (B9) 

jz=z' (BIO) 

K=ta+zj+wk', (Bll) 



where i, j, k' are unit vectors parallel to the axes of the coordinate system (x, y, z f ). Let 
the spherical polar angles in this system be 6', </>. Then 

/tan 0=tan0'. (B12) 

Now the irregularities are isotropic with respect of u, v, w, which means that all directions 
of the vector K in the system x, y, z' are equally probable. Hence, the probability p(6' , <j>)dd'd<t> 
that 6', (j), lie in given ranges dd'd<$> is given by 

p(f, 0) = (1/4tt) sin 6'. (B13) 

The corresponding probability that d,<t> lie in ranges dd,d<i> is p(0,<t>) ddd<f> where 

P(e,4>)= P (e', 0) W = 4x(cos 4 +r sin ^ )3/2 - (Bi4) 

The average values in eq (B5) to (B9) can now be found. For example 

r * p sin 5 f 27r 

sin 4 0= —, '' , n . 2Q , 3/2 dB d<t>. (B15) 

Jo 4x(cos 2 0+/ 2 sin 2 0) 3/2 Jo 

This reduces to 



cos^ 



sm^= J L(/) = (l + ^/ 2 )(l-/ 2 )- 2 -2/ 2 (l-^/ 2 )(l-/ 2 )- 5/2 log[{l + (l-/ 2 ) I/2 }//] (B16) 
Similarly it can be shown that 

riib^:=M(/)=-|/ 2 (i-/ 2 )- 2 +/ 2 (i+|/ 2 )(i-/ 2 )- 5 / 2 log[{i + (i-y^) 1 / 2 }//] (Bi7) 

Su?d=P(f) = l-L(f)—M(f) (B18) 

^Fe=Q(f)=L(f)+M(f) (B19) 
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In particular fZ(0) = l; Z(l)=8/15; M(0)=0; M(l)=2/15. 
Figure 3 gives curves showing how L(f), M(f) depend upon/. 
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Figure 3. Curves showing how the functions L(J) and M(f) 
depend on the elongation factor f. 
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